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A conversion of thermal energy into electricity is considered in the electrically polarized graphene 
stripes with zigzag edges where the heavy chiral fermion (HCF) states are formed. The stripes 
are characterized by a high electric conductance G e and by a significant Seebeck coefficient S. 
The electric current in the stripes is induced due to a non-equilibrium thermal injection of "hot" 
electrons. This thermoelectric generation process might be utilized for building of thermoelectric 
generators with an exceptionally high figure of merit ZST ~ 100 >> 1 and with an appreciable 
electric power densities ~ 1 MW/cm 2 . 

^ ■ PACS numbers: 84.60.Rb, 73.40.Gk, 73.63.Kv, 44.20.+b 

The heat energy can be immediately converted into electricity by thermoelectric generators (TEG). [H-0| However, 
available solid state TEGs typically have a low figure of merit ZST = Q oh /Q ph << 1 (here Q oh is the generated 
electric power, Q ph is the heat flow energy carried by the lattice oscillations, phonons; ST is the temperature difference 
across TEG) . [3 6] Here we consider an approach which allows improving of TEG characteristics considerably. The 
key element of the suggested thermoelectric generator (TEG) is the graphene stripe Q with atomic zigzag edges 
polarized by the transverse electric field E (see Fig. 1). The graphene TEG exploits the gate voltage-controlled heavy 
chiral fermion (HCF) states formed inside Q and it also benefits from the electric/heat current filtering. The design 
strategy is illustrated by an expression for the figure of merit ZST = G e S 2 ST 2 / A for a symmetric TEG. [3, 0] Here 
S is the Seebeck coefficient, G e is the electric conductance, A = A c + A p h is the thermal conductance due to the 
electrons (e) and phonons (ph). In typical TEG devices, |6( A p h/A c ~ 10 3 — 10 4 , providing that A ~ A p h >> A c . 
Since normally A >> G e S 2 5T 2 , one gets ZST<<1 . In the TEG reported here, instead of regular electrons and holes, 
we exploit "heavy" HCF particles to boost the electric conductance G e oc y/m*/m e and Seebeck coefficient iS (see 
Supplement). In contrast to a pristine graphene, the HCF mass m* in C/-TEG exceeds the free electron mass m e by 
a huge factor m*/m e ~ 10 2 — 10 6 >> 1. We also implement special multilayered "hot" metallic electrodes with an 
aim of decimating A p h. It allows reducing of the phonon part Q ph the heat energy flow by factor ~ 1CP 2 — 10~ 3 . The 
combined approach yields Z5T>1, or even ZST » 1, which manifests a significant improvement over existing TEG 
devices. 

The whole TEG circuit represents a 2D array formed by three intcrdigitatcd combs Q, Ga, and Gb as shown in 
Fig. Id. The comb Q consists of parallel graphene stripes alternating with metallic split-gate stripes Ga(b) which 
form the other two combs. The Ga and Gb combs serve as split-gates to creating of the transverse electric field E ^ 
across Q. Additionally, there are three other metallic combs, C , Ch and H, which are oriented in perpendicular to the 
former Q and Ga,b combs. The two combs, C e and Ch, act as heat sinks, while the H comb represents a heat source. 
Besides, there are two more back gate combs, Gl and Gr, which are placed underneath of the dielectric substrate 
and which are also oriented in-perpendicular to the Ga.b and Q combs. 

Let us consider an elementary cell of the TEG device sketched in Fig. la (we call it <7-TEG) and its diagrams 



shown in Fig. lb. The £-TEG cell consists of two field effect transistors, [8|-|l2j left (FET L ) and right (FET R ), both 
being fabricated on the same section of the graphene stripe Q. The two transistors, FETl and FETr , are connected 
electrically in a sequence, but thermally in parallel (see Fig. le). [I| The edges of each ^-stripe have an atomic zigzag 
shape. 0] The opposite edges of the same stripe are terminated by carbon atoms belonging to two distinct sublattices 
A and B. [lj| The split-gate voltage Vsg = ^ga — Vqb polarizes the C/-stripe in the transverse y-direction. The left 
(Gl) and right (Gr) back-gate electrodes are controlling the concentration and type of the charge carriers inside the 
left (x < —Lo/2) and right (x > L /2) Q-TEG sections, which are forming FET L and FET R , correspondingly. The 
local gate voltage Vq (x) induces an x-dependent change U (x) of the electrochemical potential along the (J-stripe as 
shown at the bottom of Fig. lb. 

The sharp distinctive HCF levels at energies E p = ±A result from reflections at the atomic zigzag edges 0] of the 
graphene stripes Q. The CJ-stripes are polarized in the transverse y-direction by a finite electric field E / which 
controls the value A. The origin of the HCF resonances is illustrated in terms of a mere mean field approach [il l22j]. 
The method 0, HH describes the electronic excitation spectrum E p of the graphene stripe in terms of the Dirac 
equation for bipartite sublattices. One adds a symmetry breaking term oc A into the Hamiltonian 

V. = -ihv ((a x ® i) d x + (a y <8> f z ) d y ) + U (x) (l ® i) + A (a s ® f z ) (1) 
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FIG. 1: Color online, (a) The Q-TEG cell which consists of a graphene Q stripe with zigzag edges deposited on a dielectric 
substrate. One edge of Q is terminated by atoms of A sublattice while the opposite edge by the B sublattice atoms. The 0-stripe 
is polarized in {/-direction by the transverse electric field |E| = Vsg/W = (Vga — Vgb)/W (W is the spacing between Ga and 
Gb) which is created by the split gate electrodes Ga and Gb- Two back gate electrodes Gl and Gr having opposite polarity 
induce a step- wise change of the electrochemical potential U(x) versus coordinate x. (b) The energy diagram (top panel), the 
heat and electric current diagram (central panel), and the U(x) profile along the ^-direction inside Q (bottom panel), (c) The 
multi-layered H electrode where the phonon density of states F q {ui) in adjacent layers is mismatched, (d) The TEG device 
consists of the interdigitated split-gate combs Ga(b) ° n ° ne hand and graphene comb Q on the other hand. The H comb serves 
as the heat source while the G e and Gh combs represent the heat sinks. The back gate electrode combs Gl and Gr provide 
the dynamic, coordinate-dependent n- and p-doping of Q. (e) The heat (Q) and electric (J) currents in the C/-TEG involving 
the field effect transistors (FETl.r). 



where v = 8.1 x 10 J m/s~ c/300 is the massless fermion speed, <r, and fk are the Pauli matrices, <E> is the Kronccker 
product, {i,fc} = 1...3, and A = e|E|JV/2 = eV SG /2 = e(V GA - V GB )/2, V sg is the split gate voltage, W is 
the stripe width. The pseudospin polarization [l6| means that the electric charge is depleted in one zigzag edge 
while is accumulated in the other. Then an electric dipole is formed as soon as Vsg = 2A/e ^ 0. At Vsg = 
there are two flat bands giving rise to a large density of states at the Fermi energy and being associated to zigzag 
edge statesQ- When Vsg ^ 0, the resulting pseudo-spin polarization driven by the A term in Eq. (JTJ yields an 

excitation spectrum e± = ±v\J k 2 + q 2 + (A/v) 2 ± U characterized by the energy gap 2A. The electron wavefunction 
'F satisfies the boundary conditions * A (0) = *a' (0) = at y = and * B (W) = * B ' (0) = at y = W for 
zigzag atomic edges. It gives the transversal quantization condition (QC) in the form k = q p / tan (q p W) where 

q p = ±\J {E p =F U) 2 /v 2 — k 2 — (A/v) 2 . The QC actually determines the localized energy levels when is being solved 
in respect to the excitation energy E p . Thus, for a finite A / 0, the graphene zigzag stripes Q are band insulators 
with the pseudospin polarization. The steady state energy spectrum of the C/-stripe features two fiat bottom bands 
E p = ±A which represent the highest occupied and lowest un-occupied bands. The flat bands correspond to a very 
high effective electron mass m* = (10 2 — 10 6 )m e >> m e . 

The sign of U(x) determines whether the HCF charge carriers inside of each FETl,r are electrons or holes. In 
Figs, lb, the Vg (%) profile is step-wise, so Vg (x) = —Uo/aQ for x < —L /2 (inside FETl) and Vg (x) = Uo/ac 
for x > Lq/2 (inside FETr). Here Uo is the shift of the electron electrochemical potential and «g is the back gate 
efficiency. The 5-stripe RR section located immediately under the central hot H-electrode (—Lo/2 < x < Lq/2) 
remains neutral, since V*g (x) = 0. This part of the C/-stripe serves as an electron- hole recombination region (RR). 
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FIG. 2: Color online, (a) The electron energy subbands in the graphene stripe Q with zigzag edges. The transverse split 
gate voltage Vsg induces an energy gap 2 A = eVsG = e(VciA — Vgb) which splits the narrow zero-energy level, (b) Density 
of electron states N E for three inelastic scattering rates 71,2,3 = 0.01, 0.001 and 0.0001 in units of A. The two sharp peaks 
constitute the Heavy Chiral Fermion (HCF) excitations arising at energies ±A. The corresponding peak height N e= ±a depends 
on 7. It exceeds the Van Hove singularity (VHS) peaks -/Vvhs by the big factors N e= ±A/Nvns = y/ m*/m e = 10, 100 and 1000. 
(c) The phonon transmission coefficient £ for multilayered metallic electrode H (shown in Fig. lc) versus the number of layers 
N. (d) The transmission coefficient T^ R ^ S (e) between the RR and uncovered Q sections for different electron subbands. 




Color online, (a) The transmission coefficient T(e) for 



conventional electrons and conventinal heavy fermions through the well and potential barrier, (b) T(e) for chiral fermions: 
curves 1, 2 correspond to regular chiral fermions while curves 3, 4 to the HCF excitations. Curves 1, 3 are for the incident 

angle ip = n/8 while curves 2, 4 are for ip — 3n/8. (c) The electron distribution functions in the cold electrode (C), 
recombination region (RR), and in the hot electrode (H). Inset shows the driving factor T H /rr(e). (d) The driving factor 

Yrr/c(e)- 



In the TEG shown in Fig. Id, the heat is supplied via the "hot" H-comb, while it is drained by the "cold" C c and 
Ch sink combs. The phonon fraction of the heat flow is filtered out by the multilayered H electrode (see Supplement). 
Since the electron temperature Tjj in H is much higher than the electron temperature in the (J-stripe, Tg, one achieves 
the non-equilibrium thermal injection Q ch of the "hot" electrons and holes from H into RR (see Supplement for 
details). Although the contact Seebeck coefficient <Sh/rr an d contact electric conductance Gr/rr through the H/RR 
interface vanish (i.e., 5h/rr = and Gjj/rr = 0, as shown in Supplement), the contact thermal conductance A H / RR 
remains essentially finite. [6llll| After being thermally injected from H into RR, the electrons and holes are quickly 
(during time ~ 10~ 13 s) converted into the HCF excitations. In the RR region, the non-equilibrium HCF electrons 
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and holes populate the levels ±A inversely: the upper +A level is populated by excessive HCF electrons while the 
lower —A level by the excessive HCF holes. The quasiparticle distribution function n RR is obtained as a solution of 
the quantum kinetic equation [25| (see Supplement). 

In the configuration shown in Fig. 1, the holes ballistically propagate from RR toward Ch while the electrons in 
the left (/-section proceed from RR toward C . Thus, the latter stage implies a chiral transmission of the excessive 
non-equilibrium HCF electrons from the upper +A level localized in RR into the upper A — Uq level located in the 
uncovered (/-section adjacent to C e (see the diagram at the top of Fig. lb) where Uo < 7, 7 is the localized level 
width. Simultaneously, the excessive holes are transmitted from the lower —A level in RR into the lower —A + Uo 
level localized near C/j. In this way, the full thermal flow Q eh from H to RR is eventually split between the C^ and 
C e sections of the graphene stripe Q. In the uncovered (/-section, broadening 7 of the HCF level Eq originates from 
coupling of the HCF states to the phonons. It yields 7 = 7 P h = h/r p h (typically r p h — 10~ 12 s at T = 10 K). In the 
RR section, there is an additional coupling to the electron states in H which gives 7 = 7 P h + 7h- 

Both the electrons and holes are supplied into the electrically neutral RR by the thermal injection from H. Then 
the electric current J in FETr actually emerges due to the HCF holes ballistically propagating from the RR region 
toward the C^ electrode. Simultaneously, the electric current J in FETl also consists of HCF electrons ballistically 
propagating in the opposite direction from RR toward C e . In this way, both the HCF electrons and holes transfer 
the heat from the hot (H) to the cold (C C; h) electrodes. Then the combined thermal flow Q ch of electrons and holes 
generates an electric current accompanied by a finite voltage drop Vsd = J2 n ^sb • All * ne 5-stripes are connected 
in parallel, thus the net electric current through the whole TEG is / = J^k ^ where Ik is the current through the 
k-th (/-stripe. The voltage drop on the two FETl,r is Vg D = 2S n 8T n . Because the heat flow is transferred by heavy 
charged HCF particles (m*/m e >> 1) it ensures a considerable value of Seebeck coefficient S n (see Supplement). 
It in turn yields a high thermoelectric voltage Vg^ = 2S n ST n generated by the heat flow when a finite temperature 
difference ST n is maintained between the H and C (h) electrodes of the n-th Q-TEG. Then the Q-TEG generates the 
electric power Qg_ TEG = J2k n ^^sb which is considerably higher than in TEG devices where conventional electrons 
and holes are engaged. 

Another increase of ZST is accomplished when the phonon part of the heat conductance A p h is strongly reduced. A 
visible decimation of A p h by a few orders of magnitude is achieved, e.g., if the hot H-electrodes are metallic multilayers 
fabricated as sketched in Fig. lc. The multilayercd structure shown in Fig. lc serves for decimating of the phonon 
transport through the TEG by the factor ( which we plot in Fig. 2c versus the number of layers N (see Supplement) . 
Because the multilayered metal contacts H and C e ,h are reducing the phonon component of the heat flow, the figure 
of merit ZST can additionally be increased by orders of magnitude. 

One computes the net electric power generated by the whole TEG with equivalent (/-stripes and equivalent (/-TEGs 
as Qeiectr = Ng x IVsu = Ng x G e {2S5T) 2 where Ng is the total number of Q-TEG elements in the TEG array and we 
have omitted index n. The corresponding thermal power is Qhoat = Ng x AST. Both the quantities, Qcioctr and Qhoat, 
are expressed in terms of Seebeck coefficient 5, temperature differences ST, electric G e , and thermal A conductivities 
of the individual Q-TEG sections. The TEG conversion efficiency is then estimated as {Z5T}teg = Qciectr/Qhcat = 
4:NgS 2 G e 5T/ A. The electron part, A e , of the thermal conductance A is obtained as 1/A = 2/A° + I/Aq 0011 * and is 
typically 10 3 — 10 4 times smaller than the phonon part, which is 1/A p h = 2/AP h + 1/Ag h + l/Agio 2 - [1] Here A's indices 
e, ph, c, G, and Si02 relate it to the electrons, phonons, contacts, (/-stripe, and dielectric Si02 substrate respectively. 
There are only two conducting HCF channels per each (/-stripe. For one conducting channel in a regular graphene 
stripe by width W = 10 nm, one gets Seebeck coefficient at most S ~ 10~ 4 V/K. For a graphene stripe with zigzag 
edges by the same width W ~10 nm and the split gate voltage Vsg =0.1 V, where the HCF resonances are invented, 
one improves it to S ~ 10~ 3 — 10 _1 V/K (see Supplemet). In combination with the multilayered "hot" electrode it 
yields as much as ZST ~ 10 2 . 

The electric power generated by a single Q-TEG is evaluated as Q ch ~ KHCFGgVgQ ~ 1 mW, where khcf = 
y/ m* I m e (Wh/rr/ Wrr/c) , G q = 2e 2 /h = 7.75 x 10~ 5 and we have used m*/m e = 10 2 (which corresponds to 

7 = 1 meV), <5T h /rr = 300 K, <5T RR/C = 30 K, V SG = 0.1 V. Let us assume that the length of a single £-TEG clement 
is Lq ~ 2.5 /im, the ^-stripe width is Wg = 10 nm, and the period in the y-direction is W p = 40 nm. Then the TEG 
device by size 1 cmxl cm contains Nq = 10 9 C7-TEG elements which might generate the electric power about 1 MW. 
One can see that the Q-TEG efficiency largely depends on the maximum value of Vsg and on reducing of A p h- For 
practical realizations, instead of d.c. gate voltages, one can operate the TEG at industrial a.c. frequency / = 60 or 
50 Hz. Then required Vq and Vsg might be generated immediately during the heat to electricity conversion. 

The obtained result suggests that the C/-TEG devices which exploit filtering of the electric/heat currents and of 
using the voltage-controlled spectral singularities are able to improving of the thermoelectric conversion efficiency and 
the generated output electric power by 2-3 orders of magnitude as compared to presently known devices. 
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SUPPLEMENTARY INFORMATION 

Blocking the phonon flow by multilayered electrodes. 

The goal here is to reducing of the phonon transport considerably while stimulating the electron transport at the 
same time. We design the filter pads to separate the G-section both from the external electrodes and from the 
substrate thermally, but not electrically. We consider two different types of the heat/electric current valve pads. One 
design involves metallic multi-layers Pb/Al with the layer thickness ~ 10 nm. Another method is to depositing of 
pads made of SrHfC>3 and/or SrRuOs. The layered materials have an appreciable electric conductance while their 



thermal conductance along the c-axis is remarkably low. [23|, [24| Planting of the H/RR pad between the metallic 
electrodes and graphene stripe would reduce the effective A p h significantly, because the phonons which provide large 
heat conductance A£ h between the hot and cold electrodes are eliminated from the thermoelectric path. Then the 
net heat conductance which involves the path HOT => RR => Q =>■ C c ,h => COLD will be considerably diminished. 
Thus placing of the H/RR pad with a sufficient number of nanolayers allows to decimating of the phonon part A p h 
the whole thermal conductance A. Optimal C?-TEG geometry is also determined by the electric and thermal transfer 
lengths which are estimated @, EH correspondingly as L c \ ~ 10 — 100 nm and L t h — 100 — 250 nm. 



Molecular dynamics simulations. 

The phonon part of the thermal transport through the TEG had been examined as follows. We describe the non- 
equilibrium thermal transport through the ^/-stripe in presence of multiple scattering on lattice defects, boundaries, 
and electrons. A finite temperature difference 5T between the opposite ends of each C?-stripe induces the thermal flow 
given as a sum of contributions of the individual phononic subbands. The phonon density of states Fp(ui) related to 
the phonon subband /3 is mismatched in adjacent layers of the H electrode sketched in Fig. lc. Inside the (/-stripe, 
the phonon distribution function N{uj) is non-equilibrium which means that N(lo) deviates from the Bosc-Einstein 
distribution in the hot (H) and cold (C) ends. For a "clean" graphene stripe, the phonon mean free path exceeds the 
(/-stripe length L. Therefore the non-equilibrium effect does not influence the final results. The equilibrium phonon 
distribution at the (/-stripe ends is established due to a free phonon diffusion into the bulk of attached metallic 
contacts and dielectric substrate. The thermal conductance Aq of the (/-stripe had been computed by using the 
phonon density of states -Fg(w) preliminary obtained for each phonon subband (3. 

The thermoelectric characteristics are found by solving the Dirac equation for chiral fermions in graphene (see 
above). The analytical model is verified by numeric calculations based on the density functional theory. Q The 
electron and phonon excitation spectra arc obtained considering influence of the inelastic clcctron-phonon and clastic 
electron-impurity scatterings. They are taken into account along with processes of the electron tunneling through 
the interface barriers. The electron-impurity and electron-phonon scatterings are included within the Kcldysh-Grccn 
function technique [25| which allows deriving of the quantum kinetic equations (see below for details). 



Non-equilibrium thermal injection. 

The heat energy flow Q eh from the "hot" electrode H with temperature Th into a much colder heat sinks C e n^ 
proceeds along the (/-stripe with temperature Tg (Th >> Tg). The process happens in two stages (i) and (ii) as 
illustrated by the energy diagrams at the top of Fig. lb. In the diagram we show the distribution function of "hot" 
excessive quasiparticles Sn c ± in H where 5n c = [n(e/Tn) — n(e/Tg)]8(e — £f) corresponds to the thermally excited 
excessive electrons while Sn^ = — [n(e/Tn) — Ti(e/Tg)]8(sF — e) describes the excessive hole thermal excitations [here 
n(x) = l/(e x + 1)]. (i) The thermal injection H— >RR: The temperature difference ST = Th — Tg initializes the 
non-equilibrium "thermal injection" of "hot" electrons and holes from H into the two levels ±A localized inside 
the recombination region RR. As a result, the upper +A level becomes inversely populated by the non-equilibrium 
electrons while the lower level —A is populated with non-equilibrium holes. The process is described by the quantum 



6 



kinetic equation [25] 

-q£ = ^n/RR{£,n £ } + C cp {e,n £ } + C ch {e,n e } + C RR/Ca h {£,n e } . (2) 
In the quasistationary case, one sets dn e /dt = which gives 

<nRR ^RR 

1 H/RR - n e J + 1 RR/C h (™ e - V J 0, (3) 

7"eh ^~ep 

where we have used £h/rr = Th/rr {nf — ri RR ) for the non-equilibrium thermal injection, £ cp = (n RR — n^) /r cp 
is the electron-phonon collision term, £ c h = (n RR — n^) /r e h describes the electron-hole recombination, [14] and 
£rr/C c h = r RR /c o h ( n f R — n e' h ) accounts for the HCF electron/hole escapes from RR into directions of C e ,h- 
The corresponding electron distribution functions are approximated as = 1/ (e e ^ Tli + l) in the H electrode and 
n e,h = \j ( e (e±u )/T c _|_ ^ m tne C e> h shoulders where Uq = Vo/a is the back gate-induced shift of the electron 
electrochemical potential, a is the back gate efficiency, is the Fermi function. Besides, Ih( c ) are the effective 
electron temperatures in the H and C e ,h-adjacent regions. The non-equilibrium distribution function n RR in the RR 
region then is 



_ r RR/Cc h <' h - r H / RR n" - r E n* 
r R R/c o h - r H /RR - r £ 



(4) 



where the electron energy broadening is T £ = h/r c h + h/T cp . Typically one gets [11] r H / RR = 5 meV for a rough 
G/Pd- interface while r H / RR = 60 meV for a smooth Q /Pd interface, r RR yc = 100 meV, T e = 1.5 mcV. The 
temperatures are taken as Tjj = 630 K, T RR = 330 K, and Tg = 300 K, the level position A = 50 meV, which 
corresponds to the split gate voltage Vsg = 0.1 V. The above formulas allows computing of the electric and thermal 
currents. Both types of the currents are inhomogeneous in vicinity of the H and C c ,h electrodes on the corresponding 
spatial lengths L c \ ~ 10 — 100 nm and L t h ~ 100 — 250 nm. [6,11] The "hot" conventional electrons coming from H 
in RR are converted inside RR into the "heavy" HCF excitations during the short time [11] t c ~ h/T (in Ref. [11] 
r c ~ 10 -13 s, due to the energy level broadening because the tunneling coupling between H and Q. One finds 
r = /w'hTh/rr/ (4df[) — 5 meV where 7h/rr is the H/RR interface transparency, «h is the Fermi velocity in H). 
(ii) During the next chiral transmission process, RR— > Q, which occurs on the longer timcscalc r csc < L/v ~ 10~ 12 s, 
most of the HCF electrons and holes inside RR are captured by the adjacent FETl,r. Simultaneously, minor fractions 
of HCF electrons and holes annihilate with each other [14] during the time r e h > 10 -12 s [typically r osc ~ (3 — 7) • T e h] ■ 
It sets a requirement to the spatial dimension of RR as Lq = VT e h < 1 fJ^a. Since the thermal injection is essentially 
a non-equilibrium process, the transforming of electron states between the stages (i) and (ii) is incoherent. 

Knowing n RR , one computes the electric conductance, G e , Seebeck coefficient, S, and the electronic part of the 
heat conductance, A e . The electric contact conductance Gh/rr across the H/RR interface vanishes (Gh/rr = 0) 
because the electron part of the electric current is compensated by the hole part. The H/RR contact Seebeck 
coefficient also vanishes since it is 5 h /rr = — (1/eT) (L^/L^) — Vh/rr/Wh/rr = where Vh/rr is the bias 
voltage and <5T H / RR is the temperature difference across the H/RR interface. Here we have introduced auxiliary 

functions Le = (2e 2 /h) J de ■ (e — /i) Q M e (e) T e (e) T H / RR (e) where T H /R R (e) = —d [nf — n RR ] jde is the driving 
factor, a = 1 ... 3, M e (E) = Ng (E) (Hv/L) is the number of modes, Ng (E) is the electron density of states shown 
in Fig. 2b, L is the Cf-stripc length, T e (e) is the contact transparency. Then h/ rr = 0, because the driving factor 
T H / RR (£r) vanishes as it is an even function of s. Unlikely to G r / R r and <Sh/rr, which vanish at the H/RR contact, the 
thermal contact conductance Ah/rr = (-^h/rr — [^h/rrIV^h/rr)/^ 2 *^) = L^^/ (e 2 5T) remains essentially finite. 
Using the model form M e (e) — yjm* /m e -8(e — A) one finds A h /rr, ~ J"h/rr \/m* /m e (2A 2 / h) / ST H i RR where 8T H / RR 
is an effective electron temperature difference across the i?/i?i?-interface, and -Fh/rr = (— d[nf — nf- K ]/de) e=A w 0.1. 

Along the 5-stripe, the TEG parameters 5, G e , and A e are determined purely by the electron and hole transport. 
The underlying physical mechanism is the chiral transmission of the HCF electrons and holes from the neutral RR 
section to the voltage p- and n-doped ^-sections. In the same approximation, one evaluates the electric conductance 
of FETl,r along the C/-stripe between the RR and C ej h as 

2 2 f 

Grr/c = L nn/c = -Y J rfeM ( £ ) 7RR/c c(h) 0) t rr/cO) 

-7r R /c , h , (e)T RR/c (e)| £ =A (5) 
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where TpjR/c^) is shown in Fig. 3d. Analogously, one finds Seebeck coefficient along the ^-stripe 5rr/c — 
(A/e)/ST RR/c = V SG /6T RR/G and A RR / C = L^ R/c /e 2 T - TS 2 G e ~ 2^J^T e (A 2 /5T RR/C ) /h - 2^7^ 
(A 2 /<5T RR / C ) /h—>0. The last result indicates that Seebeck coefficient <Srr/c could be huge while the electron/hole 
part of the thermal conductance A RR y c along the stripe is typically low. The phonon part of the heat energy flow is 
Q RR t c = A p h(5T RR y c = N p hKo5T RR /c where = 5 x 10~ 12 W/K at T = 10 K while the number of phonon modes 
A^ph also depends on the temperature and the stripe geometry. 
The electron/hole heat energy flow is 

n ch crp2 o2 r XT 2 o2 r _ 2e2 T /2 ^ r H/RR. V fmF 2e 2 2 

Q = 61 ■ b e{h] G e(h) - 61 ■ £>rr/ c Grr/c ~ — V SG ■ y g T c J y ~ - ~ V SG- k hcf (p) 

where we have used Vsg = 2A/e, <Srr/c = Vsg/$T rr /c, Grr/c = \/itl* /m e (2e 2 /h), and we have defined the factor 

khcf- Because khcf can be big, one might achieve huge values of Q ch . Typically A^/ t RR « A^^°, therefore the 

C7-TEG net "phonon" heat conductance [6] is A p h ~ ApY RR , which can be comparable to the contact electron heat 

conductance A H / RR = 2A 2 ^Jm* / m e FYi/ RR / {h ■ <57h/rr)- It means that only the contact electron/hole and phonon 
heat conductances actually contribute into A. Summarizing the above estimates one arrives at Z8T » 1. 

Transparency of the H/RR interface 

is directly related to the thermal injection efficiency. The interface barriers which contribute into Ar/rr originate 
from the difference of the workfunctions in the metallic H electrode and the graphene ^-stripe right beneath of it 
[11,15]. The heat-conducting Ch, H, and C e electrodes are deposited at the top of the (/-stripe, as schematically shown 
in Fig. 1. Another factor is the change in number of the conducting channels when electrons and holes tunnel from the 
3D metallic H electrode into the 2D graphene ^-stripe. [21] Conversion of the regular electrons and holes into the HCF 
excitations also contributes to Ar/rr. Thus, for the Q-TEG, Ar/rr depends on the 3D/2D conversion efficiency rj and 
on the spatial distribution of charge carriers near the H/RR interface. The contact thermal conductance problem and 
its solution are illustrated in Fig. 3. In Fig. 3a we plot the transmission probability T (e) as a function of the electron 
energy e for the conventional electrons penetrating a non-chiral potential barrier (curve 1), and the quantum well 
(curve 2). Curve 3 shows T (e) for the non-chiral heavy fermions transmitting via a potential well. One can see that 
T(e) is strongly suppressed in the latter case. For such a reason, the contact conductance for conventional "heavy" 
fermions is low. Quite a different behavior 7"(e) takes place if instead of the conventional "heavy" electrons there 
are the "heavy" chiral (HCF) particles as is evident from Fig. 3b. In Fig. 3b we compare T(e) for the conventional 
chiral fermions penetrating the chiral barrier (curves 1 and 3) with the same characteristics for HCF particles (curves 
2 and 4). One can see that T (e) is fairly good for both types of the chiral particles if the incidence angle is small, 
i.e., ip = 7r/8 (curves 1 and 2). For bigger incidence angles, i.e., <p = 3ir/8 (curves 3 and 4), for the HCF particles 
T(e) becomes suppressed (curve 4). The electron thermal conductance A e of £?-TEG is determined by 2A h /rr. The 
dominant contribution into A H / RR comes at the low angles tp, therefore using of the HCF particles helps to maintaining 
of A e at a decent level. 
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